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The multiphoton ionization of hydrogen atoms in a strong elliptically polarized microwave field 
exhibits complex features that are not observed for ionization in circular and linear polarized fields. 
Experimental data reveal high sensitivity of ionization dynamics to the small changes of the field 
polarization. The multidimensional nature of the problem makes widely used diagnostics of dynam- 
ics, such as Poincare surfaces of section, impractical. We analyze the phase space dynamics using 
finite time stability analysis rendered by the fast Lyapunov Indicators technique. The concept of 
zero-velocity surface is used to initialize the calculations and visualize the dynamics. Our analysis 
provides stability maps calculated for the initial energy at the maximum and below the saddle of 
the zero-velocity surface. We estimate qualitatively the dependence of ionization thresholds on the 
parameters of the applied field, such as polarization and scaled amplitude. 

PACS numbers: 32.80.Rm, 05.45. +b, 42.50.Hz, 92.10.Cg 
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INTRODUCTION 



The multiphoton ionization of the hydrogen atom 
in a strong microwave field is an example of a low- 
dimensional quantum system that manifests behavior 
typical of classical nonlinear systems, such as chaos and 
high sensitivity to field parameters. Ever since the first 
experiments on the ionization of hydrogen atoms in a 
strong microwave field [l| this system has been studied 
extensively for various regimes of the applied field: From 
high- to low-frequency regimes, and from linear to circu- 
lar polarization of the field [2| . The ionization dynamics 
turns out to be strongly dependent on the parameters 
of the applied field, such as frequency, field amplitude, 
the shape of the field envelope, and the choice of initial 
ensemble of states subjected to the field |3|,|5(- 



Experiments |3|, |f| on the ionization of hydrogen atom 
in strong elliptically polarized (EP) microwaves have 
demonstrated the dynamics in this case to be quite com- 
plex, exposing new effects that are absent from the cir- 
cularly polarized (CP) and linearly polarized (LP) field 
limits. The ionization dynamics was observed to be ex- 
tremely sensitive to the polarization degree of the electric 
field. Moreover, three-dimensional classical simulations 
of the experiment, reproduce the experimental results 
very well 0, • It was shown in that the ionization 
yields, deduced experimentally and numerically, interpo- 
late unevenly between LP and CP limits and for the EP 
fields follow a non-monotonic rise with the increase of 
the scaled field amplitude. By viewing the system in the 
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frame rotating with the time-dependent angular velocity 
of the EP field the authors in 0, used a quantum- 
mechanical approach to show that the problem can be 
reduced to a Rydberg atom subjected to two effective 
static fields, magnetic and electric. The sensitivity of the 
ionization to the ellipticity degree observed experimen- 
tally was ascribed to dependence of the effective magnetic 
interaction on the ellipticity degree. While the use of the 
quantum-mechanical approach for hydrogen in an EP 
field reproduces the experimental trends, a connection 
with the classical approach, which has been traditional 
in this subject, is desirable. The question still remains 
open: Why are the current experimental data on ion- 
ization yield curves of hydrogen atom in EP microwave 
fields in good agreement with the corresponding classical 
simulations |M ? 

Most classical treatments in the literature use some 
simplifications of dynamics based on the existing symme- 
tries or reduction of the full three-dimensional d yna mics 
of the system to a lower-dimensional one 

taiuiiia To 

give an example, the adiabatic approximation used by 
Griffith and Farrelly [Tl| to study three-dimensional hy- 
drogen atom in EP field reduces the dimensionality of the 
effective phase space and limits the dynamics to that of 
the circular Rydberg states. The averaging that is per- 
formed in 7] is valid only in the low-frequency regime 
(the driving frequency of the field is much smaller than 
the unperturbed Kepler frequency of the problem). The 
presence of a time-dependent term in the rotating-frame 
Hamiltonian is another challenge that the EP field brings. 
In addition, the absence of conserved quantities, such as 
energy and angular momentum integrals, prevents one 
from reducing the dimensionality of the dynamics. For 
the LP and CP problems such a reduction is possible due 
to existence of additional constants of motion. In con- 
trast, for the EP field problem the energy and angular 
momentum are not conserved. In this paper we discuss 
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the possibility to view a full-dimensional system dynam- 
ics by means of fast Lyapunov Indicator (FLI) fields of 
stability pjj . The FLI method is independent of the di- 
mensionality of the system. It provides pictures of global 
stability structures for any given subspace of the phase 
space of the system 0. Besides, the method is robust 
because it is applicable even in regimes which cannot be 
studied by any other existing analytical methods. 

The complex features of ionization dynamics of hydro- 
gen atom in EP microwave field reported in recent exper- 
iments 3] are due to the higher dimensionality of the EP 
problem as opposed to the CP and LP problems. Clas- 
sically, for systems of three or higher degrees of freedom 
the phase space dynamics is known to possess much richer 
features than the one for the lower- dimensional systems. 
The principal feature that distinguishes the phase space 
dynamics of Hamiltonian system of dimension three and 
higher is the possibility of a well-known phenomenon 
of Arnold diffusion (diffusion along the resonances) [T^ . 
This phenomenon is essential for examining all the pos- 
sible scenarios that lead to the classical ionization for 
hydrogen subjected to the electromagnetic field. 

The relevance of the phase space dynamics to the clas- 
sical ionization of hydrogen in non-stationary fields has 
been noted by many authors [1 El El EH- I n short, 
invariant tori can prevent chaotic diffusion of trajecto- 
ries and, therefore, can hinder them from ionizing. The 
onset of stochastic ionization for the Hamiltonian sys- 
tems has been traditionally investigated by means of 
the resonance overlap criterion proposed by Chirikov 
poj . It was applied to the classical ionization dynam- 
ics for the hydrogen atom driven by LP and CP field 
in @, El El HI El Hi- The set of action-angle vari- 
ables required for this analysis is only well-defined in 
the limit of vanishing field. For the classical system of 
hydrogen atom subjected to a strong EP field the appli- 
cation of action-angle variables in an overlap criterion is 
not practical. Instead, we propose to study the onset of 
stochasticity in the phase space of the system by means of 
fast Lyapunov Indicator (FLI) stability analysis [l3| . By 
evaluating the indicator of chaoticity for each integrated 
trajectory from an ensemble of trajectories we obtain the 
underlying resonant structures for any given subspace. 

The experimental data || do not provide much infor- 
mation about the character of the states that determine 
the observed ionization threshold. In fact, most of the 
initial states used in the experiment are switched by the 
turn-on of the field to different locations in the phase 
space 0, E] • Some of the states that remain bounded 
after the rise of the field pulse are closer to stochastic 
ionization than the other states. The character of the 
states more favorable to ionization has been a matter of 
debate in the theory community. Investigations on the 
ionization by the CP field in high-frequency regime were 
carried out by Howard [23| and by Sacha and Zakrzewski 
|25j . Using essentially the same technique, the Chirikov 
overlap criterion, these two groups arrived at different 
conclusions regarding to the type of the orbits that ionize 



first under the application of the field. Howard claimed 
that most orbits become elongated before they ionize, 
hence that the high eccentricity orbits are more prone to 
ionization. However, results in |25| contradict this point 
of view and the authors show that ionization occurs first 
for the medium eccentricity orbits. In Refs. 0, ^| it 
was pointed out that the ionization mechanism strongly 
depends on the region of phase space that the original 
ensemble will be switched to by the application of the 
field. Generally speaking, orbits with different eccentric- 
ities can be switched to different regions of the phase 
space and to different energies in the rotating frame of 
the CP field problem. In this paper we describe in details 
the character of the states that undergo ionization and 
determine the behavior of ionization yields for the low 
amplitudes of the field. The principal difficulty is the di- 
rect comparison of our results with experiments |3(. It 
is essential to point out, that our qualitative analysis is 
performed for ensembles of states with the same initial 
energy and different values of classical action and angular 
momentum, as opposed to the ensemble of states used in 
the experiment. Our purpose in this choice is to relate 
the underlying phase space structures to the ionization 
dynamics of the system. 

The recent experimental success in creating field- 
maintained non-spreading electronic wave packets 3 has 
renewed the discussion on synthesizing coherent states 
launched from the Lagrange equilibria of the effective 
potential (SB, It was argued in Ref. that the 

size of the stability region surrounding the Lagrange equi- 
librium points is important in maintaining stable non- 
spreading three-dimensional coherent quantum states at 
those points. The addition of a magnetic field to the CP 
problem was shown to play an important role in the sta- 
bilization of the equilibria and enlarge the critical region 
of stability in the parameter space defined by the Jacobi 
constant, and the amplitudes of the electric and magnetic 
fields. We extend these findings to the EP problem by 
using fast Lyapunov Indicator stability technique. We 
show that the stability of the Lagrange equilibria can be 
controlled by manipulating the field polarization in ad- 
dition to the Jacobi constant, and the amplitudes of the 
electric and magnetic fields. 

Much of our understanding of complex details of the 
classical ionization mechanism depends on a thorough 
knowledge of the multidimensional dynamics of the sys- 
tem. Therefore, we will apply the FLI method so as to 
obtain some insight into the stability of the phase-space 
structures of the system. 

The paper is organized in five sections. In Sec. |H] we 
introduce the Hamiltonian and describe the concept of 
time-dependent zero- velocity surface (ZVS) in the frame 
rotating with the frequency of the applied field. A brief 
description of the method of Fast Lyapunov indicators 
(FLI) is given in Sec. IIIII We discuss the results of the 
stability analysis rendered by FLI method in Sec. IIVI 
By mapping out the value of the FLI for each trajectory 
from the configuration space we obtain the FLI stabil- 



3 



ity plots. First, the structures of the FLI stability plots 
are described for the well-studied cases of LP and CP 
electric fields and zero magnetic field. The energy of the 
initial states is equal to the maximum of the ZVS. In the 
case of the CP problem the FLI plot matches perfectly 
the Poincare surface structure. Secondly, we perform the 
FLI analysis for more complicated problems of hydrogen 
in EP fields (the magnetic field is again equal to zero). 
The ionization dynamics for the EP problem is analogous 
to the LP and CP problems: two distinct sets of bounded 
orbits located around the center and near the maximum 
of the ZVS are identified. With the increase of the am- 
plitude of the electric field the orbits from the latter set 
become chaotic and ionize. The size of the stability zone 
around the ZVS maximum is used to estimate qualita- 
tively the behavior of the ionization probabilities versus 
scaled amplitude of the field. Then, the FLI stability 
analysis is applied to the ensemble of initial states with 
the energy much below the saddle of the ZVS. Section llVl 
concludes with the detailed description of two ensembles 
of initial states involved in our calculations. A summary 
of our results and conclusions are presented in Sec. 

II. HAMILTONIAN 

The Hamiltonian for a hydrogen atom subjected to 
an EP electric field (of magnitude F and microwave fre- 
quency uu) simultaneously with the static magnetic field 
B = Bz applied perpendicular to the plane of polar- 
ization is, in atomic units (ao = ft=e = /i = l) and 
assuming infinite nuclear mass, as follows: 

H = \(vl+V 2 y+P 2 z )-\B(xp y -yp x ) 



-\-F(x cos uut + ay sin uut) , (1) 

where r = \J x 2 + y 2 + z 2 and a G [0, 1] is a polariza- 
tion of the electric field : a = 1 for the circularly po- 
larized field, and a = for the linearly polarized field. 
In this paper we study the planar classical model of hy- 
drogen restricted to the x — y polarization plane. Most 
of the previous classical and quantum calculations were 
performed in the planar limit |E El El El El The pla- 
nar limit of 3D system is a reliable approximation to the 
actual three-dimensional dynamics for the orbits with 
initial conditions initiated on the z = 0,p z = sub- 
space. This subspace is invariant. Moreover the classical 
phase space dynamics for the planar limit of the hydrogen 
atom inside CP and LP fields can be studied by means 
of Poincare surfaces of section, which in this case is two- 
dimensional. Since the applied electric field is polarized 
in the x — y plane the ionization of hydrogen is expected 
to occur along any direction on the plane as in the Stark 
effect. 

Insight into the dynamics can be obtained by consider- 
ing the system in the rotating frame precessing with the 



frequency of the field uu. By introducing the canonical 
transformation 

S(x,y,p x ,p y ) = -(x cos uut - y sin uut)p x 

— (x sin uut + y cos uut)p y , (2) 

one obtains the Hamiltonian in the rotating frame (for 
convenience we omit the bars in the following expression 
and thereafter): 

H = \{pl+pl)-{l + \B){xp y -yp x ) 
r 8 

+F(a — l)(xsin 2 uut + ysinuut cos uut) . (3) 

It is more convenient to work with scaled frequencies 
and field strengths. We introduce the scaling of time, 
coordinates, momenta, and field amplitudes as follows: 
t ujt,r = uu 2 ^r,p u-^ 3 p,K = uu^^H.F' = 
uu~ 4 / 3 F,B — uu~ x B. After dropping the primes, the 
scaled Hamiltonian becomes: 

K = \{pl+P 2 y)-^ + \B){xp y -yp x ) 

-- + \b 2 (x 2 +y 2 )+Fx 
r 8 

+F(a — 1) (x sin 2 t + y sin t cos t) . (4) 

For the CP field (a = 1) the Hamiltonian is time- 
independent and the system has two effective degrees of 
freedom. The Jacobi constant introduced in celestial me- 
chanics |2?| |3(j is equal to the energy of the system in 
the rotating frame. For the EP field (0 < a < 1) the 
Hamiltonian is a time-dependent quantity. For the EP 
problem due to the absence of integrals of motion the sys- 
tem cannot be reduced to the two-dimensional one as was 
done for the CP problem. Instead, the time-dependent 
Hamiltonian in (4) can be made autonomous by in- 
troducing canonical transformation S{x 1 y 1 p Xl p yi t 1 w) = 
wK{x 1 y 1 p Xl p yi t) where w = t is a generalized coordi- 
nate. The corresponding generalized momentum is de- 
fined as p w = —K(t). The transformation yields the 
expression for the effective autonomous Hamiltonian: 

H eff = ^(pl+P 2 y)+Pw-(l + l,B)(xp y -yp x ) 

-i + ^V + y 2 ) 

+F(a — I) (x sin 2 w + y cos w sin w) . (5) 

Equation (4) describes an autonomous Hamiltonian sys- 
tem with three degrees of freedom . In practice, the 
Hamiltonian is regularized by changing coordinates to 
the set of semi-parabolic coordinates. This substitution 
is used in order to avoid the singularity near the origin 
(x,y) = (0,0) M 

An additional complexity in the problem comes from 
the the velocity-dependent Coriolis term proportional to 
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xp y — yp x in the expression for Hamiltonian (5). In its 
presence, the Hamiltonian can not be split into a pos- 
itive definite kinetic term depending on momenta alone 
and potential energy term depending exclusively on posi- 
tions. Because of the Coriolis term, the potential energy 
surface cannot be used for understanding the stability of 
equilibria. Instead, the concept of zero velocity surface, 
which constitutes an effective potential, is adapted from 
celestial mechanics [2<| |3(j . To define the Z VS we express 
Hamiltonian (4) in terms of velocities and positions : 



K 



1/.2 , -2N (1 + 

2 {X + y) ~ — 



B) 



(x 2 + y 2 ) h Fx 



+F(a — 1) (x sin 2 t + y cos t sin t) 



(6) 




Setting velocities in the above expression to zero we ar- 
rive at the following form of the effective time-dependent 
potential surface: 

V(x,y,t) = -±-<±±*l(a? + y*)+Fx 
r 2 

-\-F(a — l)(x sin 2 t + y cost sin t) . (7) 

For the CP case (a = 1) the ZVS is time-independent in 
the rotating frame. There are two equilibria lying on the 
x-axis (y = 0) corresponding to the maximum (+) and 
saddle point (— ) of the surface: 

-(1 + B)x\ + F[l + (a - 1) sin 2 t)x 2 ± ±1 = 0. (8) 

The expression for the maximum K max and the saddle 
K sa( i of the ZVS results: 



-^max(sad) i i %± 

\x±\ 2 
+F(l + (a- l)sm 2 t)x± . (9) 

Fig. [I] shows the zero- velocity contour plot and the loca- 
tion of the equilibria x±(0) at time t = 0. 

Unlike the maxima of the potential surface that are 
always unstable, the maxima of the ZVS need not to be. 
In fact, linear stability analysis carried out in the vicinity 
of equilibria in Ref. [2^ demonstrates that in the CP 
limit the stability of equilibria of the ZVS depends on 
the parameters of the field such as a scaled amplitudes 
of the electric and magnetic field. The equilibria for the 
CP problem play the same role as the Lagrange points 
L4 and L5 for the restricted three-body problem [29.l3(]|. 

Equation (7) shows that the ZVS of the EP problem 
is not constant and oscillates around the origin with the 
frequency of the applied field. The linear stability anal- 
ysis of equilibria performed in 28] for CP problem is not 
tractable for the EP problem. We study local dynamics 
around the Lagrange maximum by applying linear sta- 
bility analysis rendered by the fast Lyapunov indicator 
(FLI) method. 



FIG. 1: (Color online) Zero velocity surface contour plot at 
time t = for the hydrogen in the EP microwave electric 
field. F — 0.117, B — (in scaled units). The circle (cross) 
indicates the location of the maximum (saddle) point of the 
ZVS. 



III. DESCRIPTION OF THE FAST LYAPUNOV 
INDICATOR METHOD 

To characterize the multidimensional phase space 
structures a method is needed that provides a clear repre- 
sentation of the chaotic and regular regions in the phase 
space, similarly to the pictures obtained by Poincare sec- 
tions. However for the three- and higher-dimensional 
systems the construction and visualization of Poincare 
sections might be very complicated. An alternative is to 
define an indicator of chaoticity for each trajectory from 
a given subspace of phase space. Various diagnostics for 
different types of dynamical behavior are present in the 
literature; for instance, indicators of chaoticity are the 
characteristic Lyapunov exponent 31], the strength of 
diffusion of the instantaneous frequency in the frequency 
space 32] , the measure of complexity of the Fourier spec- 
trum j33| , and other available indicators that distinguish 
chaotic and regular dynamics 0, ^ • As a short-term 
dynamical diagnostic a method of Fast Lyapunov indi- 
cators (FLI) was proposed in Ref. [l3| to study weak 
chaos in high-dimensional systems. The main idea of 
the method is to study the evolution of the tangent vec- 
tors computed along a given trajectory of the system. 
This method has been applied to search for phase space 
structures in different types of problems: coupled stan- 
dard maps |3^, celestial mechanics problems 37], and 
vibrational dynamics of polyatomic molecules 14]. In 
this paper we apply this method to the atomic system. 
Below we present a brief description of the method (for 
more details, see Ref. |36|). 

To obtain a numerical estimate on the growth of tan- 
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FIG. 2: FLI curves for chaotic and regular trajectories of the 
hydrogen atom in a CP microwave field. F = 0.117, Kmax = 
— 1.3807, B = 0. Initial coordinates (x,y) for the correspond- 
ing trajectories are: a) (1.5, 0.01), b) (0.9,0.01), and c) (0.8, 
0.01). The values of the initial momenta are chosen as de- 
scribed in Sec. II V ATI The time t is scaled and dimensionless. 



gent vectors along the flow the equations of motion are 
integrated together with the corresponding equations for 
Jacobian matrix, i.e., we integrate the flow given by f(X) 
and the tangent flow: 



d_ 

dt 



X = f(X) , 



dj_ 

dt 



= Vf(X)J 



(10) 



where X = (xi, #2, • • • •> %n) is a n-dimensional vec- 
tor in the phase space, J is n x n Jacobian matrix 
of the flow, Vf(X) is the matrix of partial derivatives 
{dfi/dxj}ij=i^ n of the flow, and J7(0) is the identity 
matrix. The columns of Jacobian matrix are the tangent 
vectors of the flow {vj}j = i )U . Given initial tangent vector 
vo for each initial condition xo from the phase space of 
the system the fast Lyapunov indicator (FLI) is defined 
at time t as follows: 



t; x ) = Q m £ a ^ log ||v(t'; x ) 



(11) 



where v(t / ;xo) is the tangent vector at time t' . For any 
initial choice of the tangent vector vq it will eventually 
converge to the direction of the most unstable eigenvector 
of the Jacobian matrix J . Hence, characteristic dynam- 
ics can be observed by integrating Eqs. (10) for one of the 
tangent vectors from the set of all the tangent vectors in 
J. 

The Lyapunov indicator defined in Eq. (11) serves as 
an efficient diagnostic of chaoticity for any analyzed tra- 
jectory with initial condition xo. The FLI procedure sup- 
plies information about the dynamical behavior of a tra- 
jectory, analogously to the method of the characteristic 



Lyapunov exponents. The tangent vector growth eval- 
uated by integration of Eq. (10) is proportional to the 
growth of the distance between two neighboring trajec- 
tories in phase space. In fact, the evolution of the tangent 
vector computed for chaotic trajectories obeys an overall 
exponential law. Therefore, the FLI increases linearly in 
time. At the same time, the tangent vector growth for 
regular trajectories follows approximately a global linear 
behavior. Thus the FLI increases logarithmically in time. 
Already after a short interval of time it is possible to 
make a clear distinction between different dynamical be- 
haviors of trajectories by looking at the maximum value 
attained by the FLI. To provide an example, in Fig. 
the FLI evolution in time is shown for two chaotic and 
one regular trajectories from the phase space of the sys- 
tem defined by Hamiltonian (5) . Already after a time 
t = 100 two different types of dynamical behavior can 
be clearly identified. Meanwhile the curves (a) and (b) 
evolve almost linearly, the curve (c) follows logarithmic 
growth with time. Note that (b) and (c) are launched 
very close and hence the FLI are very close for some time 
before they are distinguished since (b) is chaotic and (c) 
is regular (by inspections of their Poincare sections) 

The main advantage of FLI technique over other avail- 
able methods is its power to discriminate regular from 
chaotic motions over a relatively short time period. 



IV. NUMERICAL RESULTS 



A. Ionization dynamics at the maximum of 
zero— velocity surface 



1. Choice of initial conditions 

We analyze the ionization dynamics of hydrogen atom 
in EP field by applying FLI stability analysis for trajec- 
tories from a mesh of initial conditions from the x — y 
plane. For the CP problem we consider the subspace of 
initial conditions defined in Ref. [28( • Namely, the initial 
energy is chosen to be equal to the maximum of the ZVS: 
K(0) = K max as defined in Eq. (9). The initial momenta 
Px,Py are chosen to satisfy the relation p x x + p y y = 0. 

For the CP problem the above defined subspace of ini- 
tial conditions is two-dimensional and coincides with the 
two-dimensional Poincare section. Similarly to the CP 
problem the dynamics for the EP problem is visualized 
on the x — y plane. For the EP problem the effective 
Hamiltonian of the system has three degrees of freedom. 
In order to reduce the configuration space of the system 
to two-dimensional in five-dimensional energy subspace 
one needs to specify three initial conditions for coordi- 
nates and momenta. One condition is the same as for 
the CP problem. Additional conditions are defined by 
taking generalized coordinate w = and initial general- 
ized momentum p w = —K(0). 
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2. FLI analysis for CP and LP fields 

In this section the FLI stability results are presented 
for the LP field case. The magnetic field is zero. For 
each initial condition from the subspace defined in the 
previous section Eq. (10) is integrated for a time t = 100 
(the time unit is one picosecond). The value of <j>(t) in 
Eq. (11) was evaluated at each moment of time. The 
maximum value attained over the integration interval 
[0, 100] is mapped onto the FLI plot. In Fig. 01 FLI con- 
tour plots are shown for several distinct amplitudes of 
the scaled LP electric field. The color code is assigned 
according to the maximum value of the FLI evaluated for 
each trajectory. The dark (blue in color version) color 
corresponds to the low values of the FLI and hence regu- 
lar behavior; meanwhile the light (yellow and red) color 
indicates high values of the FLI and hence chaotic behav- 
ior. All the trajectories that ionize quickly are discarded 
and not marked on the plot (white regions). In addition, 
strongly chaotic trajectories with the value of the FLI 
greater than the critical value (j) c > 50 were discarded. In 
most cases strongly chaotic trajectories with 0(100) > 50 
escape to infinity and ionize over the finite-time interval. 

In the FLI contour plots presented in Fig.|3|two islands 
of stable motions can be clearly distinguished. In fact, 
a similar structure of the FLI plots exists for amplitudes 
of the electric field in the interval [0.117,0.2]. The cen- 
tral and right islands are located at the center and at 
the maximum of the ZVS respectively (see Fig. Q). The 
structure of the central island does not change with in- 
creasing amplitude of the field. This island is constituted 
by the bounded states that remain stable in a given range 
of the field amplitude. A layer of chaotic motions is lo- 
cated around the edges of the main stable islands. By 
comparing the FLI plots on Fig. 02 with the zero-velocity 
contour plot in Fig. [I] we observe that the stable island on 
the right hand side is located exactly in the vicinity of the 
Lagrange maximum of the ZVS. The area and the struc- 
ture of the island change with increasing amplitude. The 
island remains almost the same for the field amplitudes 
F = 0.117 and F = 0.13. It shrinks at a field amplitude 
F = 0.16 and completely disappears at F = 0.2. These 
changes are the consequence of the break-up of invariant 
tori within the island. With increasing amplitude of the 
field they open paths for chaotic orbits to escape from 
the vicinity of the Lagrange maximum and to ionize. 

In Fig. Q] the FLI plots are shown for the CP prob- 
lem. First, one observes that in all the four plots the size 
of the stability island around the Lagrange maximum is 
much smaller than the size of similar structures for the 
LP problem. Second, the size and structure of the cen- 
tral island look exactly the same as for the LP case. For 
the CP problem, the bounded states inside the central 
island remain stable for all the field amplitudes from the 
interval [0.117,0.2]. The field amplitude for which all 
the trajectories launched from the right hand side island 
ionize (F = 0.15) is lower than the corresponding am- 
plitude (F = 0.2) for the LP case. As was mentioned 




FIG. 3: (Color online) FLI contour plots for hydrogen in a 
LP microwave field. Initial energy is equal to the maximum 
of the ZVS Kmax = -1.3807, B = 0. The scaled amplitudes 
of the electric field F are: a) 0.117, b) 0.13, c) 0.16, and d) 
0.2. 
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FIG. 4: (Color online) FLI contour plots for hydrogen in a 
CP microwave field. Initial energy is equal to the maximum of 
the ZVS Kmax — —1.3807, B — 0. The color scale is assigned 
according to the values of the FLI. The scaled amplitudes of 
the electric field F are: a) 0.117, b) 0.13, c) 0.14, and d) 0.15. 




FIG. 5: (Color online) The Poincare surface of section and 
the magnification of the small island from Fig.^Ja) are shown 
on panels (a) and (b) respectively. 



before, for the CP field the system is time-independent 
in the rotating frame. The structures of the FLI contour 
plot coincide with the ones of the Poincare section. As 
an example, in Fig.[5|the Trojan bifurcation is shown on 
the Poincare section ((a) panel) and on the FLI contour 
plot ((b) panel). Figure [5^6) is the magnification of the 
small island of Fig. Ufa). A one-to-one correspondence 
of both plots is evident. On a close inspection, the sym- 
metry about the y axis is observed. The resonant zones 
on the Poincare section correspond to the dark (blue in 
color version) islands on the FLI plot. The chaotic zones 
are marked around the edges of the islands. 

To illustrate the FLI stability results, we show the ion- 
ization dynamics for the trajectories with the initial con- 
ditions inside two stable regions identified on the FLI 
plots and one trajectory with initial conditions inside the 
chaotic region. In Fig. the configurations of three tra- 
jectories are shown as well as ZVS surface in the rotating 
frame. Initial conditions are chosen from different re- 
gions in Fig. 0] (a) from the big stable island, (b) from 
the small stable island, and (c) from the chaotic region. 
Trajectories (a) and (b) correspond to the low eccentricity 
bounded states and trajectory (c) is chaotic and ionizes 
after several close encounters with the nucleus. 
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FIG. 6: (Color online) Two stable orbits (a), (6), and one 
ionizing orbit (c) versus the zero-velocity surface. The pa- 
rameters are K — K max ,F = 0.117, B = 0, and a = 1. 



3. FLI analysis for EP field 

The FLI plots in Fig. [3 illustrate the structures for 
the hydrogen atom in EP field for the polarization de- 
gree a = 0.6. The scaled amplitude of the electric field 
is taken within the interval [0.117,0.17]. In general, the 
structure of the FLI plots for the EP field appears to be 
similar to the stability plots for the CP field case. Two is- 
lands of stability are apparent: a big island located at the 
center of ZVS and a small island located at the maximum 
of ZVS (see Fig. [Q. The ellipticity of the field destroys 
the symmetry that is present for the CP field case. It 
is evident that the right hand side island in Fig. (c) 
is not symmetric around the y axis as opposed to the 
ones around the maximum in Fig. EI Lightly-colored re- 
gions can be seen on the edges of the two islands. These 
are initial conditions corresponding to chaotic trajecto- 
ries with high values of the FLI. Trajectories with the 
values of the FLI higher than the critical (j) c = 50 are 
not represented. It can be observed that the large area 
surrounding the stable islands corresponds to rapidly ion- 
izing chaotic motions. Similarly to the structure of the 
central island observed in Figs. 01 and 21 the equivalent 
structure for the EP field remains essentially unchanged 
with increasing strength of the field. At the same time 
noticeable changes are seen in the size and structure of 
the small island. The size of the small island varies sub- 
stantially with the increase of the field amplitude. The 
growth of the island at certain amplitudes is accounted 
for the stabilization of some resonant motions within the 
island. For example, at F = 0.13 the size of the island is 
noticeably larger than its size at F = 0.117 in Fig. 0(a). 
At F = 0.17 the small island disappears, i.e. most of the 
invariant tori within the island that existed at lower am- 



plitudes of the field have been broken and almost all of 
the trajectories in the vicinity of the maximum become 
chaotic and ionize. The FLI plots shows that the dynam- 
ics in the EP case is more regular than the CP case. For 
example, the size of the stable zone around the maximum 
point in Fig. [7| appears to be larger than the size of the 
stable zone in Fig. 21 On the other hand the size of the 
right hand side island in Fig. [7| for the EP field is much 
smaller than the size of the corresponding island in Fig. OH 
for the LP field. From our detailed examination of the 
FLI results for the initial energy equal to the maximum 
of ZVS we conclude that the phase space dynamics is 
more regular for the LP field case than for the EP and 
CP field cases. Moreover, the size of the stable region 
around the Lagrange maximum point is observed to vary 
non-monotonically with the increase of the field ampli- 
tude for the EP field case. These changes are ascribed 
to various nonlinear resonant effects that were studied 
classically in Ref. 7] and quantum mechanically using 
Floquet state approach in Ref. Q . 



B. Study of ionization probability curves 

It has been pointed out in Ref. [01 that the appar- 
ent ionization threshold for Rydberg atoms in the CP 
microwave field must be determined by the fraction of 
orbits that undergo first transition to chaos. Such or- 
bits are located outside of the ZVS and they represent 
the atomic states that can be easily populated in the ex- 
periment |3|. In connection with these results we found 
the set of bounded orbits in the vicinity of the Lagrange 
maximum of the ZVS. These are the orbits that undergo 
the transition from a regular to chaotic behavior within 
the range of electric field amplitude F £ [0.12, 0.18]. The 
FLI plots in Figs. 03 SI an d which are computed for 
the linear (a = 0), circular (a = 1), and intermediate 
(a = 0.6) polarizations, illustrate the changes of the sta- 
bility of these orbits. While the dynamics in the cen- 
tral island remain unaffected by the increase of the field, 
the orbits located near the Lagrange maximum become 
chaotic and ionize with the increase of the field ampli- 
tude. For this reason, we determine an apparent ioniza- 
tion threshold by measuring the fraction of chaotic or- 
bits from the phase space volume enclosing the Lagrange 
maximum. The calculations of the ionization probabili- 
ties are carried out for the polarizations of the field from 
the circular to the linear limit (0 < a < 1). The behavior 
of the ionization probabilities versus scaled electric field 
amplitude F G [0.12,0.18] is estimated by means of the 
FLI analysis. 

By monitoring the evolution of the FLI along each in- 
tegrated trajectory we count the number of chaotic tra- 
jectories with the FLI value equal or above the critical 
value (j) c = 50. All chaotic trajectories are considered as 
possible candidates for the classical ionization. Indeed, 
the absence of invariant tori outside of the two stable 
structures distinguished by the FLI analysis shows that 
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FIG. 7: (Color online) FLI stability plots for hydrogen in 
EP microwave field. Ko = Kmax,ot = 0.6, 5 = 0. The scaled 
amplitudes of the electric field F are: a) 0.117, b) 0.13, c) 
0.15, and d) 0.17. 



all chaotic trajectories escape from the vicinity of the 
Lagrange maximum to a zone associated with the high 
classical action values. The FLI plots in Figs. 13141 and 
[3 show that the phase space surrounding the island at 
the Lagrange maximum is chaotic. There are no stable 
structures (except the big stable island at the center) that 
can prevent the escape of chaotic trajectories and their 
subsequent ionization. 

We define the percentage ionization probability from 
the ratio of the number of chaotic trajectories N c h aot to 
the total number of trajectories N total from the volume of 
phase space surrounding the maximum point as follows: 



^on = "T= X 100%. (12) 

^total 

The resulting apparent ionization probabilities versus 
the scaled field amplitude are shown in Fig. [H] Each 
curve represents the ionization probability for different 
field polarizations a. On the top panel of Fig. [HI results 
are shown for the polarizations of the field close or equal 
to the LP limit. The curve (a) indicates the percentage of 
chaotic trajectories in the vicinity of Lagrange maximum 
for the LP field. It is clearly seen that on average the 
percentage of chaotic trajectories within the phase space 
volume around the Lagrange maximum point is below 
20%. The ionization probability (a) increases slowly with 
the increase of the amplitude of the field. It is evident, 
that for the LP limit the dynamics around the Lagrange 
maximum remains predominantly regular for the electric 
field amplitude F in the interval [0.12,0.18]. A qualita- 
tive illustration of the changes in dynamics for the LP 
limit is given in Fig. [31 From this figure one observes 
that transition from regular to chaotic dynamics occurs 
at F = 0.2. On the top panel of Fig. [S| sharp changes in 
the behavior of the probabilities (c) and (d) are observed 
for the polarizations a = 0.2 and a = 0.3 respectively. 
First, the probabilities (c) and (d) rise rapidly from the 
20% and 40% for F = 0.117 to 80% for the higher field 
amplitudes. Secondly, ionization probabilities (c) and (d) 
exhibit non-monotonic rise with the increasing field am- 
plitude. Two pronounced local maxima in the ionization 
probabilities are observed for F = 0.16, a = 0.2 and for 
F = 0.14, a = 0.3. The main feature of the ionization 
probabilities seen on top panel is the sharp transition 
from almost monotonic behavior for a = and a = 0.1 
to the non-monotonic behavior for a = 0.2 and a = 0.3. 
This shows that dynamics around the Lagrange maxi- 
mum strongly depends on the field polarization. 

On the central panel the ionization probabilities are 
shown for intermediate polarizations. The behavior of 
the ionization probabilities is non-monotonic. For in- 
stance, there is a local maximum at F = 0.13 and local 
minimum at F = 0.162 for the curve (a). The percent- 
age of chaotic trajectories in the vicinity of the Lagrange 
maximum is equal to 80%. In fact, for the field ampli- 
tude F = 0.18 and polarization a = 0.6 the percentage 
of chaotic trajectories is almost 90%, which means that 
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most of the orbits around the Lagrange maximum ion- 
ize at these parameters of the field. These results were 
shown previously on the FLI stability plots in Fig. 
Close inspection of the top and central panels shows that 
the probabilities for the higher polarizations are shifted 
to the left with respect to the probabilities for the lower 
polarizations. For example, the maximum for the (a) 
probability from the central panel happens at F = 0.13 
and the maximum for the (d) probability from the left 
panel occurs at F = 0.14. 

On the bottom panel of Fig. [SI the ionization probabil- 
ities are shown for polarizations of the field close to the 
CP limit. The fact that the ionization probabilities stay 
close to each other for polarizations a G [0.8, 1] demon- 
strates the lack of significant variations of dynamics in 
8 the vicinity of the Lagrange maximum. The probabili- 

ties change from the 80% for F = 0.12 up to 100% for 
F = 0.18. Moreover, the ionization curves follow almost 
monotonic growth with the increase of the amplitude of 
the field. For the CP limit all the trajectories around 
the Lagrange maximum become chaotic and ionize for 
F = 0.15. This is the lowest ionization bound observed. 



C. Phase space dynamics below the saddle point of 
zero— velocity surface 

In this section we discuss the FLI stability results for 
the ensemble of states with initial energy below the ZVS 
saddle point. In addition, non-zero constant magnetic 
field was applied perpendicular to the plane of polariza- 
tion. 

Previously it has been shown that the application of 
a magnetic field perpendicular to the polarization plane 
leads to stabilization of some of the phase space invari- 
ant tori for hydrogen in a CP field [l£j, UM- I n turn, the 
invariant tori create barriers to the diffusion of chaotic 
trajectories within the phase space. In this section we 
perform FLI stability analysis for the hydrogen in EP mi- 
crowave field and non-zero magnetic field. As a starting 
point, the FLI stability analysis is applied to hydrogen 
in CP field for the same parameters of the electric and 
magnetic fields. The FLI calculations are carried out on 
a grid of points from the x — y plane for magnetic field 
B = 0.2 and initial energy Kq = — 2 (the saddle point 
value of the ZVS is K sad = -1.7022). The resulting FLI 
stability plots are shown in Fig. O 

On panel (a) the structure of the phase space is shown 
for the CP field. Note the dark (blue in color version) 
island of stability around the nucleus. The same stable 
island was observed in Figs. 03 21 an d [3 It is located 
around the center of the ZVS. The island is surrounded 
by the classically inaccessible region. The phase space 
outside of the region is foliated with regular and chaotic 
motions. Chaotic dynamics is associated with high val- 
ues of the FLI and color coded with light colors (yellow 
and red in color version). The large resonant zone is lo- 
cated on the right from the center of the FLI plot. The 
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FIG. 9: (Color online) FLI stability plots for hydrogen in 
an EP microwave field. F = 0.117, B = 0.2, K = -2. The 
polarization of the field is : a) 1, b) 0.8, c) 0.5, and d) 0.3. 
The color scale is assigned according to the maximum values 
of the FLI attained over the integration interval t — 100. 



phase space appears to be mostly regular, except for thin 
chaotic layers winding around the large resonant zones. 
It is easy to notice that the structure is symmetric with 
respect to the y axis. On panel (b) the phase space dy- 
namics is pictured for the polarization a = 0.8. The 
structure resembles the one shown on panel (a): the FLI 
plot shows the same resonant structures and chaotic lay- 
ers as those observed on panel (a). However, the over- 
all dynamics appear more chaotic than for the CP field 
case. The stability results for the intermediate polariza- 
tions a = 0.5 and a = 0.3 are shown on the panels (c) 



and (d) correspondingly. The phase space dynamics on 
both panels appear to be much more chaotic than the 
dynamics for the CP limit. Although the large part of 
the phase space on contour plots (c) and (d) corresponds 
to strongly chaotic motions, several stable islands can be 
distinguished. For example, two small islands of stability 
are located symmetrically with respect to the y axis. 

Our main observation from the FLI stability results de- 
scribed in this section is that the application of magnetic 
field leads to the stabilization of the resonant structures 
within the phase space of the system. The FLI plots 
reflect the onset of stochasticity in the phase space of 
the system that occurs when the field polarization ap- 
proaches LP limit. The phase space structure appears 
to be more regular for the field polarizations close to the 
CP limit as opposed to the polarizations close to the LP 
limit. 



D. Description of an ensemble of initial states and 
comparison with the stability analysis results 

Finally, in this section we introduce the classification 
of the ensemble of orbits used for the FLI stability cal- 
culations in the previous sections. In the experiment 
the ensemble of initial states are prepared spanning a 
narrow range of the Keplerian energy E = —-^z and a 
variety of high eccentricity orbits. As the field is turned 
on, the states with the same action n will appear in dif- 
ferent regions of the phase space having diff erent values 
of energy and eccentricities e = yl — l^/n 2 (l z is angu- 
lar momentum) in the rotating frame. While the Hamil- 
tonian dynamics is traditionally studied by restricting 
the phase space to the constant energy manifold, one 
would ideally prefer to relate numerical results with the 
experiment and to study an ensemble of states with the 
same Keplerian energy. However, most of the numerical 
simulations based on the examining the structure of the 
Poincare section use trajectories with different values of 
the Keplerian energy K. For the FLI stability analysis 
we use an ensemble of initial states with equal initial en- 
ergies K(0) and different values of initial action n and 
eccentricity e. We choose these initial states to compare 
the stability pictures given by the FLI method with the 
Poincare sections. 

The interpretation of numerical results and comparison 
with experiments is much more complicated for the EP 
problem. The presence of the time-dependent term in 
the expression of the Hamiltonian in the rotating frame 
prevents the one-to-one relation between the initial state 
(at time t = 0) and the final state (obtained after the 
turn-on of the field). Therefore, it is instructive to relate 
the FLI stability results for the ensemble of states ana- 
lyzed with their initial eccentricities, actions and angular 
momenta. The importance of the configuration of initial 
states in the ionization mechanism was realized, e.g., in 
Refs. USEE^. We found that our results agree in their 
key features with the conclusion stressed out in Ref. [2]} : 
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FIG. 10: The action variable n, angular momentum l z , and 
eccentricity e of the initial states of hydrogen in the EP mi- 
crowave field are shown on panels (a), (6), and (c) correspond- 
ingly. The parameters are F = 0.117, K m ax = —1.3807, and 
B = 0. 



the fate of the initial states subjected to the application 
of the field cannot be predicted from the initial action 
values and eccentricities of such states. In essence, the 
states with the same initial action and eccentricity values 
can be moved by the application of the field to different 
parts of the phase space. We illustrate these conclusions 
by comparing the FLI stability results given in Sees. HV XI 
and II V Cl with the configuration of initial states involved 
in the FLI computations. 

First, we analyze the ensemble of states with initial 
energy equal to the maximum of the ZVS. In Fig.lTUlthe 
action variables n, angular momenta l z and eccentricities 
e are computed for the ensemble of initial states used 
for the FLI stability analysis in Sec. II V Al The case of 
zero magnetic field is considered. The scaled amplitude 
of the electric field is F = 0.117, polarization is a = 1, 
and the initial energy value is equal to the maximum 
Kmax = -1.3807 of the ZVS defined in Eq. (9). On 
panel (b) the values of initial angular momentum l z are 
shown. It can be seen from Fig. HUT M and FLI plots in 
Fig. [7| that the states inside the central stable island cor- 
respond to negative angular momentum states. The rest 




FIG. 11: The action variable n, angular momentum l z , and 
eccentricity e of the initial states of hydrogen in the EP mi- 
crowave field are shown on panels (a), (6), and (c) correspond- 
ingly. The parameters are F = 0.117, K = -2, and B = 0.2. 
The initial energy of the states K is taken below the saddle 
point energy K sa d — —1.7022. 



of the trajectories in the FLI plots in Fig. corresponds 
to positive angular momentum initial states. Another 
important observation is that the central stable island 
is associated with bounded states of different eccentrici- 
ties. From the stability results given in Sec. II V these 
states remain bounded for the range of the scaled field 
amplitude [0.117,0.2] and for different polarizations of 
the field. From Fig. HUT c) and FLI plots in Fig. [7| it is 
apparent that the initial states around the maximum of 
ZVS are the low eccentricity states. In fact, these are the 
states that determine the low ionization threshold com- 
puted in Sec. IIVBI Summarizing data given in Figs. Q 
and we argue that the ultimate fate of initially low 
eccentricity and positive angular momentum states is to 
be moved by the driving field to the chaotic region of the 
phase space and undergo fast ionization. At the same 
time, the initial states with negative value of angular 
momentum located around the center of the ZVS remain 
bounded. 

Next, we classify the initial states with energy below 
the saddle point of the ZVS. In Fig.[H]the action, angular 



13 



momentum and eccentricity together with the FLI plot 
are presented for the initial states of hydrogen subject 
to EP field F = 0.117 with polarization a = 0.9 and 
magnetic field B — 0.2. The initial energy K(0) = —2 
is chosen below the saddle point K sa d — —1.7022 of the 
ZVS. The main conclusion from the data given in Figs. El 
and is that the initial states with low eccentricity tend 
to be more stable than the states with high eccentricity 
values. For example, from Fig. HJcf) one observes three 
small stable islands with low values of the FLI: one island 
is around the center of ZVS and the other two islands are 
symmetric with respect to the ?/-axis. Two symmetric 
islands correspond to low eccentricity states in Fig. HTT c). 

The comparison of data presented in this section with 
the stability results given by FLI analysis show that, al- 
though the fate of initial states can be determined from 
their eccentricity and angular momentum values in some 
cases, in general, this information is not sufficient for pre- 
dicting the dynamics of the states after the turn-on of 
the field. Instead, the analysis of the character of the 
states should be performed in conjunction with the sta- 
bility analysis. 

V. CONCLUSIONS 

In this article, we provide a qualitative description of 
the classical phase space and its relevance to the ioniza- 
tion of hydrogen atom in a strong EP microwave field. 
Using the FLI stability analysis, the complex multidi- 
mensional dynamics is shown to depend sensitively on 
the changes of parameters including polarization, ampli- 
tude of the electric and magnetic fields and the initial 
state ensemble. The FLI stability analysis allows us to 
picture complex phase space structures, such as impor- 
tant resonant and chaotic zones. 

We map out the FLI values for each trajectory origi- 
nating on the zero-velocity subspace for the ensemble of 
orbits with initial energy at the maximum of the ZVS. 
Our FLI stability computations for the LP and CP fields 
reveal two main resonant structures: a small stable island 
around the Lagrange maximum of the ZVS and a large 
stable island around the center of the ZVS. The small sta- 
ble island corresponds to the initial states with low eccen- 
tricity values and positive angular momenta. In contrast, 
the initial states from large stable island are the negative 
angular momentum states with different eccentricity val- 
ues. These states remain bounded and their dynamics 
is almost unaffected by the application of the microwave 
field. The remaining part of phase space is attributed to 
chaotic motions that ionize quickly over the interval of 
time considered in the computations. The main feature 
of FLI stability results is that the ionization dynamics is 
determined by the orbits with low eccentricity and pos- 



itive momentum states located around the maximum of 
the ZVS. These states are the first to become chaotic and 
ionize for a given strength of the field. Similar observa- 
tions hold for the EP field case: The FLI plots reflect 
stable structures similar to those observed for the LP and 
CP problems. The main distinction from the LP and CP 
field results is the behavior of stable island around the 
Lagrange maximum. It changes non-monotonically with 
the increase of the amplitude of the electric field due to 
break-up and stabilization of resonant torus structures 
within the island. 

Our main conclusion is that the nonlinear stability of 
the Lagrange maximum is determined by the parameters 
of the field, such as the amplitude of electric and mag- 
netic field as well as the field polarization. The ioniza- 
tion probabilities versus the amplitude of the electric field 
reflect the changes in the local dynamics around the La- 
grange maximum. They manifest non-monotonic growth 
for intermediate ellipticities of the field. This behavior 
agrees with the experimentally observed ionization yields 
inRef. Q. 

The FLI stability analysis was also carried out for the 
ensemble of states with initial energy below the saddle 
value of the ZVS. In this case the FLI plots reflect the 
onset of stochasticity in phase space that takes place as 
the field polarization changes from the circular to the lin- 
ear limit. The FLI simulations reveal the detailed struc- 
ture of phase space that is more regular for the field po- 
larizations close to the circular limit as opposed to the 
polarizations close to the linear limit. 

The application of the method can be potentially im- 
portant for the control of ionization of Rydberg states 
in recent experiments |3^. The main advantage of this 
technique over other available analytical and numerical 
methods is its independence from the dimensionality of 
the system. Thus, the stability analysis can be carried 
out for any subspace of initial conditions in the full- 
dimensional phase space of the system. 

The method can be useful for estimating the size of 
a stable zone around the equilibria for the wave packet 
calculations and stabilization of quantum wave packets 
at the Lagrange equilibria points discussed in Ref. [28^ . 
We expect that by varying polarization and amplitude of 
electric and magnetic fields, one can achieve the stabiliza- 
tion of a set of resonances around the Lagrange maximum 
that can possibly lead to the stabilization of the quantum 
wave packets in the EP field system. 
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